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The massive Dirac equation on a Kerr-Newman background may be solved by the method of 
separation of variables. The radial and angular equations are coupled via an angular eigenvalue, 
which is determined from the Chandrasekhar-Page (CP) equation. Obtaining accurate angular 
eigenvalues is a key step in studying scattering, absorption and emission of the fermionic field. 

Here we introduce a new method for finding solutions of the CP equation. First, we introduce a 
novel representation for the spin-half spherical harmonics. Next, we decompose the angular solutions 
of the CP equation (the mass-dependent spin-half spheroidal harmonics) in the spherical basis. The 
method yields a three-term recurrence relation which may be solved numerically via continued- 
fraction methods, or perturbatively to obtain a series expansion for the eigenvalues. In the case 
p = zboj (where u and ^ are the frequency and mass of the fermion) we obtain eigenvalues and 
eigenfunctions in closed form. We study the eigenvalue spectrum, and the zeros of the maximally 
co-rotating mode. 

We compare our results with previous studies, and uncover and correct some errors in the litera- 
ture. We provide series expansions, tables of eigenvalues and numerical fits across a wide parameter 
range, and present plots of a selection of eigenfunctions. It is hoped this study will be a useful 
resource for all researchers interested in the Dirac equation on a rotating black hole background. 

PACS nurabers: 04.70.-s, 04.62.-|-v 



I. INTRODUCTION 

The interaction of a fermionic (Dirac) field with a rotating charged (Kerr-Newman) black hole is of deep theoretical 
interest. The interaction depends on three couplings: coupling between quantum-mechanical spin and black hole 
(BH) rotation; coupling between field mass and BH mass; and coupling between field charge and BH charge. Desire 
for a deeper understanding of these couplings has motivated a range of studies over the last three decades. 

Over thirty years ago, Chandrasekhar [T] and Page |Hj showed that separation of variables is possible for the massive 
Dirac field on the Kerr-Newman spacetime. In other words, the partial differential equations (PDEs) governing the 
evolution of the Dirac field may be reduced to a set of coupled ordinary differential equations (ODEs), and the solution 
can be expressed as a sum over modes and integral over frequency. The radial and angular ODEs are coupled through 
an angular eigenvalue A. The eigenvalue is found by solving the so-called Chandrasekhar-Page (CP) equation (see 
Eq. (|4| below) . The eigenfunctions of the CP equation - the so-called mass dependent spin-half spheroidal harmonics 
(MDSHs) - are required for a full reconstruction of the field. In this paper we present a new method for computing 
both the eigenvalues and eigenfunctions. Our aim is to show that finding MDSHs is no more difficult than finding 
the spin- weighted spheroidal harmonics for massless fields [3]. Along the way, we review alternative methods and 
highlight some inaccuracies in the literature. 

Motivation for this work arose from two separate studies, conducted independently by the present authors. The 
Kerr-Newman solution takes a very beautiful form in the limit where the electromagnetic field is assumed to dominate 
over gravity, which is obtained by letting G ^ in the metric. The result is fiat-space, in oblate-spheroidal coordinates. 



'Electronic address: 
t Electronic address: 



sam.dolan@ucd.Ie] 
j gairSast . cam . ac . uk 



2 



with an electromagnetic field 

E + iB = -V ( ^ ) (1) 

V V (i* - • (r - la) j 

i.e., a Coulomb field centred at the imaginary point r = ia = i(0,0, a) [H [SJ [3. This field can be used to model a 
rapidly rotating but effectively massless nucleus [7 and in his thesis jH] JG examined the possibility that the hyperfine 
splitting observed in muonium could be reproduced by the coupling between electron spin and the frame dragging 
caused by nuclear rotation in this model [7]. SD studied the spectrum of fermionic quasi-bound states in the vicinity 
of a small non-extremal black hole with a < M, extending recent work in this area [SJ [TD]. Both studies required 
the determination of accurate eigenvalues for the angular eigen-equation, and both studies were impeded by a lack of 
reliable numerical results in the literature. These studies will be presented in full elsewhere. 

The massive Dirac-Kerr-Newman system has been studied in many other contexts. For example, authors have 
investigated the absence of fermionic superradiance [TTl [T5]; scattering and absorption [T31 EJ 1121 HSJ [T7]; and 
Hawking radiation emission [THl [H]. Our intention here is to provide a solid foundation for possible further work. 
Let us give three examples of future applications. First, it has been shown [20l EH [22l [231 El] that a rotating BH is 
stable to fermionic perturbations. Nevertheless, the quantative effect of field mass [5S] on the fermionic quasi-normal 
mode spectrum [22] has not been studied. Secondly, field mass is usually neglected when the Hawking emission 
process is considered |2]. The emission spectrum will change significantly if the Hawking temperature approaches 
the mass of a particle species of the standard model. Emission of massive scalars by a Kerr-Newman black hole was 
considered in |27j , but a corresponding study for massive fermions is still lacking. Finally, in recent years there has 
been much interest in theories with extra dimensions which lead to the possibility of BH creation in hadron colliders 
[28] . The radial equations for fermionic fields on the brane are closely related to their 4D counterparts [29], and 
the angular equations are unchanged. The effect of field mass on emission of fermions from non-rotating holes was 
recently considered |30j . The method outlined here will be of use to study the emission of massive fermions from 
rapidly-rotating higher-dimensional black holes. 

This paper is organised as follows. In Sec. |ll| we introduce the basic equations. In Sec. |III| we examine the non- 
rotating limit, and introduce a novel form for the spin-half spherical harmonics. In Sec. |IV| we make explicit the 
symmetries of the eigenspectrum. In Sec.|V]we derive exact expressions for the eigenvalues and eigenfunctions for the 



special cases = iaw. In Sec. VI we define a spectral decomposition method for tackling the general case, which leads 
to a three-term recurrence relation. In Sec. [Vlll we use the recurrence relation to investigate the asymptotic behaviour 
of the eigenvalue spectrum in the slow-rotation regime. In Sec. |VIII| we describe some alternative methods used to 
find eigenvalue solutions, and highlight some errors in previous work. Numerical results are presented in Sec. |IX[ 
We examine the dependence of the eigenvalues on a/i and aui, investigate the zeros of the eigenvalue spectrum, and 
present a gallery of eigenfunctions. We conclude with a brief discussion in Sec. |X| Tables of eigenvalues, plus simple 
but accurate numerical fits to the data, are given in Appendix |B] 

II. BASICS 

In 1976 Chandrasekhar [1] demonstrated how to separate variables for the massive Dirac field on the Kerr spacetime 
and shortly afterwards. Page [2] extended the analysis to the Kerr-Newman spacetime. More recently, extensions to 
non-asymptotically flat [31] and other related spacetimes [31] have been considered. In such analyses, the t and 
(j) dependence may be factorised by using the ansatz 'if{t,r,6,(j)) = e~*"*"'"™'^^'ijm('', and the Newman-Penrose 
method is applied [33] . The four components of the wavefunction 'i'uim can be expressed as products of two radial 
and two angular functions {R^{r), R_{r), S+{d), S^{9)} that obey coupled ODEs. The angular equations are: 

+ ( ^ cot 6* — aojsin^ + TOCSc6' ) 5'_|_ = — (A — a^cos 6*) 5_, (2) 
du \2 J 

dS. (\ 



• 1 „ cot -I- aw sin 61 — m CSC 6* S'_ = + (A -f a^cos0) 5*+. (3) 
du \ 2 ) 

Here, a — J /M is the black hole angular momentum parameter, and lo and /i are the frequency and mass of the 
state under consideration (N.B. the symbol a = —oj is also used for the frequency in the literature [31]). S+{9) and 
S- (9) are known as mass-dependent spheroidal harmonics of spin one-half, and A is the eigenvalue. The eigenvalue 
can be regarded as the square root of a generalised squared total angular momentum [35]. The solutions depend on 
two continuous parameters, au and a^. For a given auj, a/i, the eigenstates {S+{9),S-{9),X\ may be labelled by 
three discrete numbers: the angular momentum j = 1/2, 3/2, . . ., the azimuthal component of the angular momentum 
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j and the parity V = ±1. When this dependence is to be made explicit, we will write S± =s=±i/2 
jmV ■ j^MLiauiuiiD and hence the eigenvalue spectrum and eigenfunctions, exhibit a number 

;ie spectrum in the quadrant aw > 0, a/i > is sufficient to determine the 



m = - j, + 

and A = 

of symmetries (see Sec. |IVp. Knowledge 
full spectrum. 

The first-order equations ^ and ^ may be combined to obtain a second-order equation, 



ions 

'wledge oitl: 



'Cot6'- 



i cos 6*) 2 



d0 



+ [(a^w^ - a^M^) cos^ 6 + aw cos 0] S- + afi sin 9 S-^ 



[X — a LO +2auj7ri\S- 



(4) 



known as the Chandrasekhar-Page angular equation. In the limit /i = this reduces to the well-known equation for 
the spin- weighted spheroidal harmonics (see for example Eq. (2.1) in [3 , making the identifications s = —1/2 and 
sAim = -I- 2aujm — a^w^). We find, in the case a/i ^ 0, it is actually easier to analyse the coupled first-order 
equations directly. 

There have been a number of studies of Eqs. (|2}|4|) over the years [13j [34 



I EH EH Eg. In Sec. VIII we discuss 



alternative approaches in some detail, to validate our method, and to identify some errors in the literature. 



III. EXACT SOLUTIONS FOR a = 



In the non-rotating limit, a — 0, the equations reduce to 

dS+ (\ 



d0 
dS- 



cot 6 + m CSC 9 j 
cot 9 — m CSC 9 1 S- 



-XS- 
-XS+. 



Making the substitution 



leads to coupled equations 



's+ 




■ cos(6'/2) 


sin(6l/2)' 




'T+ 


S- 




-sin(0/2) 


cos(6'/2) 







d 


'T+ 


+ 


d9 


T_ 





{m+l/2)cot9 (to+A+1/2)' 
(m-l/2-A) -(to- 1/2) cote 





'T+ 







These can be written in second-order form. 



sm 9 d9 \ d9 J sm 9 
which is simply the general Legendre equation. The eigenvalues are 



A 



(0,0) 
'jrriP 



= V{j + 1/2) 



(5) 
(6) 

(7) 

(8) 

(9) 
(10) 



where V — ±1 and j is a half-integer. The eigenfunctions are T±{9) oc P™*^^^(cos 0), where p™^^/^ [g an associated 
Legendre polynomial obeying the Condon-Shortley phase convention, and L = j + 'P/2isa, non- negative integer. Let 
us define a set of solutions 



+ l/2^'m-p(e) 
-l/2Yjmv{&) 



1 (j - m)] 
2tt {j + m)! 



1/2 r 



cos(6'/2) sin(6'/2)' 
-sin(6l/2) cos(6'/2) 



P 



1-1-1/2 



T,rfi— 1/2/ 



(11) 



where x — cos 9, Cj„iv = 'P{L + 1/2) — m and L = j + V /2. 

It is straightforward to verify that these solutions arc normalised so that 



27r / d9sii-i9,Yj^v{0)sYj^^v{0) 
Je=o 



(12) 
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and that they exhibit the symmetries 



sY,^v{0) = (13) 

= 7'(-i)'"-i/2-,y„-p(0), (14) 

= V{~iy+"' -sY,rnv{T^-e). (15) 

Note that by combining these symmetries, we may obtain a further four expressions of similar form. 

The solutions (jlll are closely related to the spin-weighted spherical harmonics of spin- weight half, s=±i/2y\m^\s), 
first introduced by Newman and Penrose [40] . We make the identification 

1/2 



' 2j + l(j + m)!(j-m)! 
Att {j + s)\ {j-s)\ 



(16) 



[sin(6l/2)]2j' 
(-l)^-''-"[cot(6l/2)]2"+''-™ 



(17) 



The representation (17 1 was introduced by Goldberg et al. (see |4T], eq. 3.1). The representation of s3^,-m^'' terms 



s (Eqs. ( [TT| ) and (16)) does not seem to be well-known in the literature. These functions arise 

and we find that the representation (111 has many 



of Legendre polynomia 

in the guise of spherical monogenics in geometric algebra 
advantages. 



IV. SYMMETRIES 



The eigenvalue spectrum has the following symmetries 

Hence knowledge of the spectrum in the quadrant aw > 0, a/i > is sufficient to determine the full spectrum. In 
close correspondence with Eq. (13)-(14l, the eigenfunctions have the following symmetries 



- 7'(-l)™-i/V.4l™^(0), 
= 7'(-l)J'+"-,5j™'(^-^), 



(19) 
(20) 
(21) 



and combinations thereof. For example, combining (19l and (20l yields sS'j^"^'(6') = 'p(^—iy^+"^ __,Sj _°^'p'^^'(^) 



V. EXACT SOLUTIONS FOR au = ±afi 
In the special case that auj — afi we find that transformation ([7| yields the coupled equations 

= 

i±l/2 



d 


'T+ 




de 




+ 



(m -I- 1/2) cot 6* {\ + m + l/2 - auj) 
(to- 1/2- A-acj) -(m - 1/2) cot 6* 



(22) 



Again, the solutions are associated Legendre polynomials, T± oc Pj^ , where L is a non-negative integer. The 
normalised solutions are 



A 



cos{e/2) sin{e/2) 
-sm{e/2) cos(6l/2) 



where L = j + V /2 and bjmv — [{L + 1 /2)^ — ni^] /(A + to+1/2 — aui) and the normalisation factor is 



A- 



_2L+1 (L- TO - 1/2)! / 



V{m — auj) 



47r {L + m+ 1/2)1 \ {L + 1/2^ - 2mauj + a'^uj^ 



(23) 



(24) 
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The corresponding eigenvalue is 



jm'P 2 



(25) 



In the special case m = — j, P = —1 we find b\rn\,m,-i = 1 and X\m\,m,-i = ~ j ~ old- 

The solution for auj — —ajj, can be found by considering the symmetry of equations ^ and (|3| under the simulta- 
neous transformations 



a/i ^ -a^l, A ^ -A, ,y,™p(0) ^ (-i)^-i/^yj,„_p(0) 



(26) 



Hence the solutions are 



auj , — auj 



auj. — auj 



jm'P 



(0) 
(0) 



B 



cos(e'/2) -sin(6l/2)' 
sin(6l/2) cos(6i/2) 



i(auj. — au>) p7n— 1/2/ 



(27) 



where now L = j- V /2 and feJCp^""^ = -[{L + 1/2)2 _ ^2]/(_^ 
factor is 



jmP "'^^ + TO + 1/2 — acj) and the normalisation 



^2 _ 2L+ 1 (L-TO - 1/2) 



1 - 



V{m — au) 



Att (i + TO+l/2)!\ ^(L+ 1/2)2 -2TOacj + a2w2 



(28) 



The corresponding eigenvalue is 



[a^^au) ^ 1 + 1/2)2 _2maa; + a2^2. 



jm'P 2 



(29) 



VI. A SPECTRAL DECOMPOSITION METHOD 



An obvious next step is to seek an expansion for the spheroidal harmonics in the basis of spherical harmonics (111. 
This was the approach followed by Chakrabarti [37] , who applied the method to the second-order Chandrasekhar-Page 
equation Q to find a five-term recurrence relation. Here we apply the method directly to the first-order equations 
themselves, to recover a three-term recurrence relation. The key advantage of a three-term relation is that it can be 
solved using robust continued- fraction methods. 

Let us begin with the ansatz 



(30) 



k' — \m\ 



where sC^.^^^p'^^^-j are expansion coefficients to be determined, and s^fc'm'p(^) were defined in ( |ll| ). The symmetric 
(151 and (21) imply that the coefficients are related by 



s''fc'm-p(j) 



(-1) 



For clarity, we neglect the least-relevant indices, writing ±i/2C^°top(!) = c^^A-, and ±i/2Yk'mV = ^t^^ ■ 



ansatz (30) into equation (^2h and using (31 1 leads to 



(31) 

Substituting 



^ (-V{k' + l/2)yir^ - sin OY^^^ + (-l)^'-^' (A - a^i cos 9)Y^7^) c[+l^ = 0. 



k' — \m\ 



(32) 



If we now multiply by Y^ and integrate we obtain the matrix eigenvalue equation 



Akk'hk' = Xbk 

A;' — |m| 



(33) 
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where 



h = (-l)^-'=cW =4^, (35) 



and 



Cill = 2tt I d0sm9Y^->{e) cos9Yy{9), (36) 



\ 1/2 

-1 (/c'lmO|fcm) (fc'l iO|fci) , (37) 



^(1) 

^kk' 



/ 2k' + 1 ^ 
\2k + l 

= 2tt d9 sin eYl~\d) sine Y^;^\e), (38) 

"'0 



/ 9k' -\-^\ ^1"^ 

= V\2 (A:'lmO|fcm) (fc'l - il|/fci). (39) 

The Clcbsch-Gordan coefficients are only non-zero for fc'— f<fc<fc' + f. Hence A].y is a tridiagonal matrix, and 
the eigenvalues A and vectors 6^ can be found via linear algebra routines. Alternatively, we can obtain a three-term 
recurrence relation for the expansion coefficients c^.^?^ . Using expressions for the Clebsch-Gordan coefficients given in 
Appendix]^ it is straightforward to show that 

afe6fc+i -I- 7fc6fe_i = 0, fc = |m|, |m| -M, . . . (40) 

where 

, V(fc-f 1)2 -m2 , ^ 

Uk = [an + ekauj) 2{k + I) ' 

, s / aiom \ aum , , , 



— 

7fc = [a^i-ekai^) — . (43) 

ZK 

and efc = (— l)-'"'^?-'. Note that j, fc and to are half-integers. The advantage of a three-term recurrence relation is 
that it can be solved via continued-fraction methods. To obtain numerical results, we followed the approach outlined 
in 133], and used a rescaling algorithm described in |35]. The basic idea is to separate out the eigenvalue by writing 
= /?fe — A, and write down an expression for the ratio of consecutive terms 

= ^Ji (44) 

^'-^ /3fc-A + «fc('^) 

The series should converge and this condition gives the possible values for the eigenvalues. In practice, we write the 
eigenvalue as a continued fraction 

^ = f^^-7> ^ "'^'a37. (45) 

This fraction can be continued to a certain level, and then we use the fact that the ratio hk+i/hk ^ as fc ^ oo to 
ignore subsequent terms. This is an equation of the form A = g(X) which can be solved iteratively, by starting with a 
guess for A, evaluating the right had side of (45 1, using this as a new estimate of A and repeating. 



VII. ASYMPTOTICS 



Below we show that the new three-term recurrence relation (40 1 may be used to study the eigenspectrum in the 
small-aw, a/x limit. To understand the opposite limit, auj,a^ — > oo, we recall some key results in the literature. 
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A. Small aco and afi 

For small values of acu and a/i one may express the separation constant A as a power series. We start with the 
continued fraction equation in the form 

^ "j-l7j / «|m|7|m|+l \ ^ Qj7j + 1 aj + ilj+2 ^^g^ 

^ Pj-l - ■ ■■ \ P\m\ J + 1 - ■ ■ ■ fij+2 - • • ■ 

and expand the separation constant as a Taylor series, 

C30 oo 

W = EE (47) 

p=0 q=0 

By grouping together like powers of auj and a/i we obtain expansion coefficients, 

Aoo = V{] + l/2) (48) 

Aio = -\VmK+{j) (49) 

Aoi = \mK-{j) (50) 

A20 - V[H{j)+H{j + l)] (51) 

A02 = A20 (52) 

A30 - \Vm[Kt{j)H{j) + Kt{j + l)H{j + l)\ (53) 

Ao3 = im[i^r(j)iJ(j)-i^r(j + 1)^^0 + 1)] (54) 

An = 2[i/(i + l)-i/(j)] (55) 

A21 = \m[{2K+{] + \)-K^{j + \))H{3 + \)-{2Kt{])-K^{3))H{])\ (56) 

A12 = imT' [(X+(j + 1) - 2ifr(j + 1)) H{j + 1) + (i^+(j) - 2K^{j)) H{j)] (57) 

where (for compactness) we have defined the functions 

H{k) = (fc2 - m2)/(8fc3) (58) 
K^{k) = l//c±l/(fc+l) (59) 
Kt{k) = l/[(fc+l)(fc-l)]±l/fc2. (60) 

If desired, the series may be continued to higher orders with the aid of a symbolic algebra package. It is straightfor- 
ward to confirm that the expansion coefficients are consistent with the exact eigenvalues for aoj = ±afi presented in 
Sec. |v]and expansions for the massless case {a^ = 0) given in, e.g., In Sec. |VIIl| we validate against an alternative 
series expansion in powers of {auj ± a/i) obtained in |34l I36j . 

B. Large auj, massless a/i = 

In the massless case (a/i = 0), asymptotic results for the eigenvalue in the large-ao; limit were obtained by Breuer, 
Ryan and Waller [46 . In our notation, 

= 2{q - m)aw + A + Axl(auj) + Aijiaujf + Aj,l{auSf + ©((aw)"'') (61) 

where 

^0 = -i {q^ -rr? + 2s + \\ (62) 

Ax ^ ~\{q^ -Tr?q + q- s^{q + m)\ (63) 

^2 = (5g^ - (6m2 - 10)(?^ + - 2n? - '^s^{(^ - - 1) + l) (64) 

and s = —1/2 and .43 is given in |46J. Breuer tt al. leave the parameter q undetermined. Casals and Ottewill j47] 
showed that q may be obtained by counting the number of zeros of the solution. We obtain the best match between 
( 61 1 and our numerical results when we take q = j + z, where z = if | j — m| is even and z = 1 otherwise (see Fig. |2| . 



However, we are unable to find a suitable asymptote for the m = j mode. 
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FIG. 1: Solutions to the asymptotic problem. Lowest modes of the solution to the asymptotic recurrence relation, Eqs ( 65 1-( 68 1 , 
for m — 1/2 (left panel) and m = —1/2 (right panel), as a function of the fixed mass-to- frequency ratio r = /i/oj. Modes with 
A(r) > have parity "P = 1 and those with A(r) < have parity V = —1. 



C. Large auj and afi 

The solutions with afi = are a special case of the asymptotic behaviour. This is most easily seen from the 
Chandrasekhar-Page equation, Eq. (|4|. When afi = 0, the terms quadratic in auj can be combined into a single term 
-a^w^sin^ 9S-. As a result, if the eigensolutions are confined to a region near 6* = of size A9 < l/y^auj, this term 



is actually 0{auj) and hence we expect ~ auj. It is clear from the results of Breuer et al. [46', given in Eq. (61 1 
that this is indeed the case for the asymptotic a/i = solutions. 

If instead we consider the limit auj — )■ oo, with the ratio r ^ fi/uj ^ fixed, there is an additional quadratic 
term, —a^uj^r^cos^9S-, in Eq. It is not possible for both this term and the —a^ uj"^ sm"^ 9S- term to be small 
simultaneously, and so we expect that A^ o^lu'^ in this case. This is borne out by our numerical results. Writing 
A = aujX{r) -I- ■ • • , the asymptotic solutions A(r) may be found by taking the limit aw oo of the recurrence relation. 



Eqs (40H43I 



akh+i + Pkbk + %h-i = 0, fc = |m|, |m| -f 1, . . . (65) 

where 

~ £fc(fc + l/2)m m ~ 

^ — Hk+i) + mrrj 

7. = (r-e^)^^^^. (68) 



and efe = {—iy~^V as before. We can solve this recurrence using the same method outlined in Section VI The results 
are shown in Fig. [ijfor < ?' < 1. We note that as r ^ 1, all the positive parity modes converge at A = 1 and 
the negative parity modes converge at A = — 1, which is consistent with the aw — > oo limit of the exact solution for 



r = 1, Eq. (25). We have verified numerically that the eigenvalues of the full problem do indeed have a linear slope 
asymptotically, and that this slope is correctly predicted by the asymptotic solutions shown in Fig. [T] 



VIII. ALTERNATIVE METHODS 



As discussed briefly in Sec.|TT] alternative methods for finding eigensolutions exist in the literature, although some 
mistakes are also present. Here we seek to review and clarify the present situation. 
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SufFern, Fackerell and Cosgrove (hereafter SFC) fM' showed that the eigenfunctions could be expressed as an infinite 
series of hypergeometric 2F1 functions, and derived a three-term recurrence relation for the series coefficients ( ^34j . Eq. 
(16)). They showed that the series terminates for the special case aoj — ~afi, and derived an exact result compatible 
with p9| ) in this case. Their three-term relation can be solved via the continued-fraction method to determine A. We 
have verified that the eigenvalues calculated via SFC's recurrence relation agree with the eigenvalues calculated via 



( 33 1 . A key result presented in SFC was a series expansion of the form 

\ = ^^Crs{aa - anYiaa + afiy. (69) 

where a = —u>. Unfortunately, the expansion coefficients presented in Table I of [34' are not all correct (for example, 
diagonal elements C„„ should be zero [36 ). In the second and fourth tables of Table II of [SJ, the eigenvalues are 
only accurate to one decimal place for m — ±1/2 and j = 3/2 and j — 1/2. 

A five-term recurrence relation was calculated by Chakrabarti |37j : however, we believe the result should be treated 
with caution because the numerical eigenvalues presented are incorrect. For example, in Tables I, Ila, lib and III in 
[37], the eigenvalues are wrong, except for the special case a/x = 0. There seems to be a tacit assumption made that 
the spectrum is symmetric about zero (i.e. that if A is an eigenvalue then so is —A), which is not the case for 7^ 0. 



The correct symmetries of the spectrum are given in Eq. ( 18 1 



Kalnins and Miller T3] also obtained a three-term recurrence relation, and a series expansion in powers of a, viz. 
A — '^r=o^'r'^^ ■ Ranganathan [48 presented exact solutions for the special case uj — ±fj,, m — ±j, which are 
compatible with those given in Sec. Iv] and in SFC [34] . 

Batic, Schmid and Winklmeier [56]ihereafter BSW) recently showed that the eigenvalues of ^ satisfy a first order 
partial differential equation, 

d\ d\ 

(p. - 2i)X) -^ + (i>- 2jl\) — + 2mjl + 2fw = 0, (70) 
dpi ov 

where jl — afi and v = —aui. Using this PDE, BSW found a method for calculating the series coefficients Crs in 



SFC's expansion (69 1. BSW point out that the method applied by SFC to determine the series coefficients is plagued 
by divide-by-zero problems at certain orders {r^s). BSW describe a method for avoiding the problem, allowing the 
series expansion to be taken to arbitrary order. We have checked that the expansion coefficients given in [36] are 



in full agreement with our expansion, Eq. ([48|-(57l. The non-linear PDE (70) may be solved by the method of 
characteristics, and BSW show that the characteristic equations can be reduced to a Painleve III equation 36J. 

Recently, Batic and Nowakowski [35] (hereafter BN) derived an ordinary differential equation for the eigenvalues 
at fixed mass-to-frequency ratio r ~ fi/cu , using analytic perturbation theory. This equation was 

dX 2(auj — mj)(2Xau} — afi) 



da a(4A2 - 1) 



(71) 



Unfortunately, this equation is not correct (see discussion below). In the corrected calculation, the eigenvalue is given 
once again by a partial differential equation, so a solution is not so easily obtained. 

We have found that it is possible to adapt the approach of BN in order to derive an alternative recurrence relation 
for the angular eigenvalues. Following the sam e notation as BN, we write 5*4. = gi / ^/sin 6 , S'_ — 32/Vsin 6, k — — m, 
rUe = fJ, and M = a. In deriving the ODE (71 1, BN introduced functions U = + §2, V — g2 — gl and W = 2gi g2, 
which obey the differential equations 

U'(e) ^ -2f{e)V(9) + 2MmeCoseW(e), f{e) ^ Mujsme+ (72) 

sm^ 

V'{9) = -2f{9)U{0) + 2XWie) (73) 

W'{e) = 2MmeCoseUi9) -2XVi9) (74) 



These are equations (3.5)-(3.7) of BN, but we have corrected a factor of 2 in equation (72 1. BN state correctly that 
dX/dM = TOg/i -I- UJI2, where 

h= cos6'y(6')d6', l2= sm0W{e)de. (75) 
Jo Jo 

BN went on to derive three equations for Ii, I2 and a third integral, I3 — cos^ 6U (6) d9 . However, the sign of uj 
on the left hand side of their third equation (number 3.15) was wrong, with the consequence that they appeared to 
obtain three independent equations, whereas in reality the third equation was a linear combination of the first two. 
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The fact that their original solution was not correct was indicated by the consequence that 1^ = 0, while I3 is the 
integral of the square of a real function and must therefore be greater than zero. 

It is possible to adapt and extend their approach as follows. We introduce two families of integrals, which generalise 
the integrals used in BN 



J n 



cos"6' 5^(6») d0, 



sm o cos 



2n , 



-gmUO) d9, 



(76) 



Due to the symmetry of the eigenfunctions 51 (tt — 9) = 52(^)1 and the fact that they are normalised such that 
gfde = 1, we see that Jq = 2, Ji = h, Li = h and J2 = h- Multiplying Eq. ^ by sixiO cos^""! 6* and 
integrating we find 



(2n - 1) J2„_2 - 2nJ2n = 2MmeLn + 2MujJ2n+l - 2{MlU + fc) J2„-l. 

Multiplying Eq. ( 73 ) by sin 9 cos^" 6 and integrating gives 

2nJ2n-l - (2n + 1) J2n+1 = 2MujJ2n+2 " 2{Muj + k)J2„ + 2XL„. 



Finally, multiplying Eq. (74 1 by cos^"+^ 9 and integrating gives 

(2n + 1)L„ = 2MmeJ2 



n+2 



2\J- 



2n+l- 



(77) 
(78) 
(79) 



These three relations together can be solved recursively to obtain all of the integrals. For eigenvalues of the system, 
J„, L„ — > as n ^ cxD, and this can be used to solve the recurrences. Combining the three equations gives a five-term 
recurrence for A 



(4^2 - 1) J2„-2 + 2{2n + 1){Mlu + /c) J2„_i - 2n{2n + 1) J2n 
+ (2(2n + 1)A - 2Muj).J2n+i ~ Ud^mlJ2n+2 = 



or a four-term recurrence that now depends on 



2n{2n + 1) Jsn-i + 2(2n + l)(Mw + k)J2n - ((2n + 1)^ - 4A2) J2„+i 

-(2(2n +1)Mlu + 4MmeA) J2„+2 = 0. 



(80) 



(81) 



It is clear from the above that it is possible to use the approach of BN to derive the eigenvalues, without resorting to 
a spectral decomposition. However, the result is a recurrence relation and the nice analytic solution for the eigenvalue 
derived in the original version of their paper is lost. While it may be possible to modify their approach further to 
obtain the eigenvalue exactly, we have so far been unable to do this. The recurrences in Eqs. (80) -(81 1 are more 
difficult to work with than the three-term recurrence in Eq. (40 1 and so we have used expression (40 1 to compute all 



the results presented in Section IX 



IX. RESULTS 



In this section we present some numerical results for the eigenvalues and eigenfunctions of equations ^ and (|3|. 
These were computed using the three-term recurrence relation, Eqs. (40)-(43l, and the series expansion in spin-half 
spherical harmonics (30 1. We have cross-checked against results from the other techniques (Sec. Villi. As before we 
separate variables using the ansatz '^{t,r,9,(j)) = e~^'^*e^"^'^^ ^ra{r,9). 



A. Eigenvalues 



Figure [2] shows eigenvalues for the massless case, = 0, for the modes j = 1/2, . . . , 11/2 and parity V ~ +1. The 
eigenvalues are split on azimuthal number m. For positive auj the lowest-lying eigenvalue is the m = j mode, which 
tends to zero as au 00. For negative au, the symmetries (18 1 mean that the spectrum looks the same but the 
ordering of eigenvalues in m is reversed (i.e. the m = — j mode is lowest-lying). It also follows from (18 1 that, in the 
massless case /i = 0, the eigenvalues of negative parity V — —1 are found by inversion, A —>■ —A. 

The lower plots in Fig. [2] compare numerically-determined eigenvalues against asymptotic results of Sec. |VII| In 
the small-aw regime, we use expansion (47 1, taken to third order. In the large-aw regime, we compare against the 
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Large aoj asymptotics, j = 5/2 
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m = +5/2 









3.5 4 



FIG. 2: Massless Spheroidal Eigenvalues (ajj, — 0). The top plot shows the eigenvalues A for modes j = 1/2... 11/2, m = — j . . . j 
in the parameter range aco = 0...8 and afi = 0. For positive aui, the m — +j eigenvalue is smallest; it tends to zero as 
au} — > 00. The massless spectrum (a/x = 0) is symmetric under A — > —A. The lower plots compares numerical results [solid] 
with approximations [dotted] valid in the small-aci; (left) and large-aci; limits, given by Eq. (47 1 and Eq. (61 1. 



asymptotics (61 1. Note that there does not seem to be a satisfactory asymptote at large-aw for the eigenvalue that 
approaches zero (e.g. m — j for ato > 0). 

Figure [3] compares the eigenvalue spectrum for the case auj — afi (Eq. ( pS] )) with the massless (a/i = 0) spectrum. 
The j — 5/2 modes are shown here; other modes follow a similar pattern. Observe that the m — j , P = —1 eigenvalue 
crosses the A = line at high auj, but no other eigenvalue changes sign. This can be understood by examining the 
exact result, Eq. (25 1. For V = —1, m = —j we obtain A = + 1/2) + aui and hence A = when au — j + 1/2. 
Likewise, the m = +j eigenvalue is zero when au = —{j + 1/2). For all other states, the term under the square root 
in (25) is non-zero for all au > 0, and hence the eigenvalue does not change sign. 

The zeros of the eigenvalue A are worth some further consideration. Figure |4] shows that, for each j, the m = j, 
P — —1 mode has zero eigenvalue along a line in the quadrant au > 0, a/i > 0. We find no evidence to suggest that 
the eigenvalue of any other m, V mode passes through zero in this quadrant. The zeros in the other quadrants follow 
immediately from the symmetries of the spectrum (Eq. 18 1. 

In Appendix |b] we present some tables of numerically-determined eigenvalues, for the j = 1/2 and j — 3/2 modes 
in the range 0<aa;<l,0</i/w<l. Our intention is to provide a resource for checking and validating future 
studies. 
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FIG. 3: Eigenvalues for j = 5/2, m = ~j ■■■j and V — ±1. The solid lines show the /i = uj eigenvalues (Eq. 25 1 for positive 
(red) and negative (blue) parities V. The dotted black lines show the massless eigenvalues jj. — 0, also shown in Fig. [2] The 



eigenvalue of the afi — au, 



5/2, V = —1 mode crosses the axis at au — 3 (i.e. A 



(3,3) 
5/2,5/2, 



_i=0). 



B. Eigenfunctions 



Figures [5] and [6] show the eigenfunctions S+{9) (red) and S-{9) (blue) of the j = 1/2 and j = 3/2 modes, for a 
range of aw > and a/x > 0. Figure |5] shows the massless spheroidal harmonics {afj, = 0). The solid line represents 
the spherical (a = 0) harmonics, and the broken lines represent the spheroidal harmonics with auj ~ 1,2,3 and 4. 
As the rotational coupling increases, the eigenfunctions show a tendency to increase in magnitude near the poles 
(6* ~ 0°, ~ 180°), and decrease in magnitude around the equatorial plane {9 — 90°). 

Figure [g] shows j — 1/2, 3 — 3/2 eigenfunctions for the special case w = /i, in the range atj = . . . 4. In this case, 
the eigenfunctions are known in closed form ( [27| . The plot makes it clear that the eigenfunctions of the V = +1, 
m = ±j, afi = auj mode do not depend on auj. For the other modes, we again see a general 'enhancement' towards 
the poles. 

The eigenfunctions plotted in Figs IS] and [6] exhibit the symmetries stated in Eq. ( 19 )-j[20|. For example, the 
reflection symmetry S+{9) = 7'(— 1)^"™S'_(V — 9) is obvious. In the massless case (Figj^there is also a clear 
symmetry under parity change V = —V; whereas in the general case, this only holds if we simultaneously flip the 
sign of the mass ~fJ-)- Solutions for arbitrary auj, afj, can be found by applying the symmetries (19H(21| to the 

solutions in the first quadrant. 

Figure [7] shows the influence of fleld mass on the shape of the eigenfunctions. Here, we plot the eigenfunctions 
of the j = 1/2, j = 3/2 modes, for mass-to-frequency ratios n/uj from zero to two. In general, it would seem that 
mass seems to lead to additional structure in the eigensolutions around the equatorial plane. However, we defer any 
attempt at physical interpretation to a future in-depth study. 



X. CONCLUSION 



We have described a technique for the determination of eigenvalues and eigenfunctions of the angular equation for 
the massive Dirac field on a rotating black-hole background. By carrying out a spectral decomposition of the angular 
eigenfunctions, we obtained a three-term recurrence relation. Exact eigenvalues can be easily and quickly obtained by 
numerically solving the three-term relation (40 1; alternatively, they may be estimated from series expansion (47 1. We 



have described two alternative methods to obtain the eigenfunctions, including an alternative set of recurrence relations 
that were derived in an attempt to correct one of the previous attempts to obtain a solution for the eigenvalues 
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FIG. 4: Contours in the acj-afi plane along which A = 0. The three lines show maximally co- and counter-rotating modes with 
j = 1/2, j — 3/2 and j — 5/2 (red, blue, green). The plot illustrates the symmetries of the eigenvalue spectrum, given in 



Eq. (181 



We have tabulated eigenvalues for a range of values of the spin of the central black hole and the mass of the fermion 
field (Tables |l - VI I , and derived fits which reproduce the numerical results to high precision (Tables VII XII ) . 

The motivation for this work was the need for angular eigenvalues as input for two separate studies — an inves- 
tigation of bound states in the massless Kerr-Newman background [8J and a study of the spectrum of quasi-bound 
states in the vicinity of a small black hole jUITO]. The solution to the radial equation for these two problems will be 
presented elsewhere, together with an estimate of the error arising from using fitting function Eq. (Bll in place of 
exact angular eigenvalues. 

We hope this work will also be of benefit to other studies, for example, investigations into (i) the effect of field mass 
on fermionic quasi-normal ringing of black holes; (ii) the effect of particle mass on the Hawking radiation emission 
spectrum for temperatures close to the mass of the particle species; and (iii) the interaction of brane-localised fermions 
with rotating higher-dimensional black holes. 
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FIG. 5; Angular solutions S+{6) [red] and S-{9) [blue] for ajj, = 0; the spin-half spheroidal harmonics. The top four plots 
show the J = 1/2 modes, and the bottom eight plots show the j = 3/2 modes, as a function of polar angle 6. The co-rotating 
modes m > are on the left, and the counter-rotating modes m < are on the right. Plots come in pairs: one for positive 
parity (upper) and negative parity (lower), V — ±1. Note the symmetry under 9 ^ n — 9. The solid line shows the solution 
for auj — 0. The broken lines show the solutions for auj — 1.0, 2.0, 3.0 and 4.0. 
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APPENDIX A: CLEBSCH-GORDAN COEFFICIENTS 



In this section we give explicit forms for the Clebsch-Gordan coefficients that appear in integrals ( 37 ) and ( 39 ) , and 



present simplified expressions for C^J^,' and I?^],' . The relevant coefficients are 



and 



and 



{k — 1, 1 m I A: m) ~ 
(fc, 1 m I A: m) — 

{k + 1, lmO\km) = - 

{k-1, liO|fci) : 
(A:,liO|fci) : 

(fc+ 1, 1 iO|fci) : 

(fc-l,l-il|fci) 
(fc,l -il|fci) 



{2k){2k~l)' 
m 



/2[(fc+l)2-TO2] 

(2fc + 3)(2fc + 2)^ 



I2k + 


1 


k 


1 1 


k{k^ 


-1) 


I2k 


+ 1 



fc + 1 



(fc + 1, 1 - illfci) = 



2fc + 1 
8fc 

2fc + 1 

y^8k{k + 1) 

I 2k + l 



8{k + l) 



Hence the integrals ( 37 ) and ( 39 1 are 



Vfc2 -TO2/(2fc) 

m/[2fc(fc+ 1)] 



fc' = fc - 1 
k' ^k 

V(fc+l)2-m2/[2(fc + l)] fc' = fc + l 
|fc-fc'|>l. 



Vfc2 -m2/(2fc) 

-m (fc + l/2)/[fc(fc + 1)] 

-v/(fc + l)2 -m2/[2(A: + 1)] 





k' ^k-l 
k' = k 

fc' = fc + 1 
l/c-fc'l > 1. 



(Al) 
(A2) 

(A3) 

(A4) 
(A5) 
(A6) 

(A7) 
(A8) 

(A9) 



(AlO) 



(All) 



APPENDIX B: TABLES OF EIGENVALUES 



Numerical eigenvalues for the j = 1/2 and j = /2 modes are tabulated in Tables |l- VI The parameter range is the 
same as in |37] . aa — Q . . .1 and r — [iju = ... 1. The eigenvalues are accurate to six decimal places, and have been 
checked using the method of SFC [34 . They disagree with the eigenvalues presented in Chakrabarti [3T, Tables I, 
Ila, lib and III. It is possible to find simple polynomial fits that reproduce these eigenvalues with high precision. We 
take an ansatz of the form, and fit it in the range auj e [—1,0], r = fi/uj e [0, 1] 



=4fe=4 



A - Ao + ^ ^ /„fe(ac^)"(l + ^i|oJf. 

n=l k=0 

The coefficients fnk that reproduce the data in Tables |l VI are tabulated in Tables VII XII 



(Bl) 
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p = 


+1 






aw 






R = 




0.1 


0.2 


0.3 


0.4 


0.5 


0.0 




0.934097 


0.869818 


0.807306 


0.746712 


0.688186 






1.067385 


1.136116 


1.206065 


1.277108 


1.349131 


0.2 




0.941098 


0.884521 


0.830461 


0.779115 


0.730674 






1.061036 


1.124021 


1.188780 


1.255145 


1.322958 


0.4 




0.948159 


0.899463 


0.854152 


0.812462 


0.774612 






1.054745 


1.112158 


1.172009 


1.234081 


1.298167 


0.6 




0.955278 


0.914640 


0.878364 


0.846712 


0.819912 








1. lUUOZo 


1. 100 / 00 




1 .Z / 4 ( ( 


0.8 




0.962457 


0.930048 


0.903081 


0.881822 


0.866479 






1.042342 


1.089133 


1.140022 


1.194663 


1.252724 


1.0 




0.969694 


0.945683 


0.928286 


0.917745 


0.914214 






1.036229 


1.077973 


1.124808 


1.176305 


1.232051 






0.6 


0.7 


0.8 


0.9 


1.0 


0.0 




0.631876 


0.577928 


0.526476 


0.477646 


0.431544 






1.422021 


1.495674 


1.569992 


1.644879 


1.720244 


0.2 




0.685319 


0.643219 


0.604522 


0.569349 


0.537786 






1.392071 


1.462344 


1.533646 


1.605854 


1.678854 


0.4 




0.740800 


0.711195 


0.685925 


0.665076 


0.648682 






1.364073 


1.431621 


1.500643 


1.570983 


1.642500 


0.6 




0.798151 


0.781565 


0.770229 


0.764154 


0.763286 






1.338022 


1.403483 


1.470931 


1.540169 


1.611021 


0.8 




0.857194 


0.854030 


0.856969 


0.865911 


0.880683 






1.313893 


1.377878 


1.444410 


1.513243 


1.584153 


1.0 




0.917745 


0.928286 


0.945683 


0.969694 


1.000000 






1.291647 


1.354724 


1.420937 


1.489975 


1.561553 



TABLE I: Eigenvalues A for j = 1/2, m = ±1/2 and P = +1. In each row, the top line is the m = +1/2 eigenvalue and the 
lower line is the m = —1/2 value. Here, R is the ratio between mass and frequency u!. 



p = -l 






aui 






R = fj,/u! 


0.1 


0.2 


0.3 


0.4 


0.5 


0.0 


-0.934097 


-0.869818 


-0.807306 


-0.746712 


-0.688186 




-1.067385 


-1.136116 


-1.206065 


-1.277108 


-1.349131 


0.2 


-0.927156 


-0.855357 


-0.784701 


-0.715290 


-0.647230 




-1.073792 


-1.148441 


-1.223859 


-1.299963 


-1.376672 


0.4 


-0.920276 


-0.841142 


-0.762657 


-0.684881 


-0.607877 




-1.080258 


-1.160994 


-1.242157 


-1.323697 


-1.405563 


0.6 


-0.913456 


-0.827176 


-0.741186 


-0.655515 


-0.570190 




-1.086781 


-1.173773 


-1.260952 


-1.348294 


-1.435774 


0.8 


-0.906697 


-0.813461 


-0.720298 


-0.627215 


-0.534219 




-1.093362 


-1.186776 


-1.280236 


-1.373736 


-1.467268 


1.0 


-0.900000 


-0.800000 


-0.700000 


-0.600000 


-0.500000 




-1.100000 


-1.200000 


-1.300000 


-1.400000 


-1.500000 




0.6 


0.7 


0.8 


0.9 


1.0 


0.0 


-0.631876 


-0.577928 


-0.526476 


-0.477646 


-0.431544 




-1.422021 


-1.495674 


-1.569992 


-1.644879 


-1.720244 


0.2 


-0.580627 


-0.515589 


-0.452221 


-0.390625 


-0.330897 




-1.453910 


-1.531603 


-1.609681 


-1.688073 


-1.766714 


0.4 


-0.531711 


-0.456447 


-0.382152 


-0.308894 


-0.236738 




-1.487709 


-1.570089 


-1.652658 


-1.735371 


-1.818185 


0.6 


-0.485243 


-0.400704 


-0.316603 


-0.232973 


-0.149846 




-1.523370 


-1.611058 


-1.698816 


-1.786622 


-1.874453 


0.8 


-0.441318 


-0.348520 


-0.255831 


-0.163261 


-0.070817 




-1.560827 


-1.654406 


-1.747998 


-1.841598 


-1.935199 


1.0 


-0.400000 


-0.300000 


-0.200000 


-0.100000 


0.000000 




-1.600000 


-1.700000 


-1.800000 


-1.900000 


-2.000000 



TABLE II: Eigenvalues A for j = 1/2, m = ±1/2 and P = —1. In each row, the top line is the m = +1/2 eigenvalue and the 
lower line is the m = —1/2 value. Here, R is the ratio between mass /u and frequency u!. 



p = 


+1 






au) 








R = 




0.1 


0.2 


0.3 


0.4 


0.5 


0.0 




1.920329 


1.841350 


1.763119 


1.685695 


1 


609143 






2.080312 


2.161215 


2.242663 


2.324614 


2 


407027 


0.2 




1.924473 


1.849942 


1.776486 


1.704190 


1 


633144 






2.076449 


2.153751 


2.231842 


2.310661 


2 


390152 


0.4 




1.928643 


1.858641 


1.790098 


1.723128 


1 


657851 






2.072612 


2.146386 


2.221237 


2.297085 


O 

Z 


o / OoOD 


0.6 




1.932840 


1.867446 


1.803952 


1.742503 


1 


683251 






z.UDo / yy 


9 1 QQ1 1 Q 


91 nS/1 Q 


Z.ZOOOOO 


2 


358139 


0.8 




1.937063 


1.876356 


1.818048 


1.762311 


1 


709329 






2.065012 


2.131951 


2.200679 


2.271067 


2 


342997 


1.0 




1.941311 


1.885372 


1.832381 


1.782542 


1 


736068 






2.061250 


2.124881 


2.190725 


2.258623 


2 


328427 






0.6 


0.7 


0.8 


0.9 


1 





0.0 




1.533529 


1.458925 


1.385410 


1.313063 


1 


241970 






2.489862 


2.573085 


2.656660 


2.740555 


2 


824737 


0.2 




1.563445 


1.495195 


1.428503 


1.363483 


1 


300256 






2.470262 


2.550942 


2.632146 


2.713831 


2 


795954 


0.4 




1.594395 


1.532895 


1.473492 


1.416335 


1 


361576 






2.451480 


2.529894 


2.609036 


2.688852 


2 


769289 


0.6 




1.626354 


1.571979 


1.520298 


1.471488 


1 


425727 






2.433514 


2.509932 


2.587316 


2.665597 


2 


744709 


0.8 




1.659291 


1.612391 


1.568825 


1.528789 


1 


492475 






2.416357 


2.491044 


2.566964 


2.644029 


2 


722158 


1.0 




1.693171 


1.654066 


1.618962 


1.588061 


1 


561553 






2.400000 


2.473214 


2.547950 


2.624100 


2 


701562 



TABLE III: Eigenvalues for j = 3/2, m = ±3/2 and P = +1. 



p = -1 






auj 






R = h/lo 


0.1 


0.2 


0.3 


0.4 


0.5 


0.0 


-1.920329 


-1.841350 


-1.763119 


-1.685695 


-1.609143 




-2.080312 


-2.161215 


-2.242663 


-2.324614 


-2.407027 


0.2 


-1.916210 


-1.832864 


-1.749998 


-1.667649 


-1.585858 




-2.084199 


-2.168777 


-2.253701 


-2.338945 


-2.424480 


0.4 


-1.912118 


-1.824486 


-1.737125 


-1.650054 


-1.563297 




-2.088112 


-2.176436 


-2.264955 


-2.353651 


-2.442509 


0.6 


-1.908053 


-1.816216 


-1.724500 


-1.632913 


-1.541465 




-2.092050 


-2.184194 


-2.276423 


-2.368731 


-2.461110 


0.8 


-1.904013 


-1.808054 


-1.712125 


-1.616228 


-1.520366 




-2.096012 


-2.192048 


-2.288106 


-2.384182 


-2.480276 


1.0 


-1.900000 


-1.800000 


-1.700000 


-1.600000 


-1.500000 




-2.100000 


-2.200000 


-2.300000 


-2.400000 


-2.500000 




0.6 


0.7 


0.8 


0.9 


1.0 


0.0 


-1.533529 


-1.458925 


-1.385410 


-1.313063 


-1.241970 




-2.489862 


-2.573085 


-2.656660 


-2.740555 


-2.824737 


0.2 


-1.504668 


-1.424125 


-1.344277 


-1.265176 


-1.186878 




-2.510280 


-2.596321 


-2.682579 


-2.769030 


-2.855652 


0.4 


-1.476877 


-1.390819 


-1.305151 


-1.219898 


-1.135092 




-2.531512 


-2.620646 


-2.709896 


-2.799247 


-2.888686 


0.6 


-1.450166 


-1.359027 


-1.268059 


-1.177273 


-1.086681 




-2.553551 


-2.646048 


-2.738593 


-2.831179 


-2.923800 


0.8 


-1.424541 


-1.328755 


-1.233011 


-1.137312 


-1.041660 




-2.576385 


-2.672507 


-2.768640 


-2.864783 


-2.960933 


1.0 


-1.400000 


-1.300000 


-1.200000 


-1.100000 


-1.000000 




-2.600000 


-2.700000 


-2.800000 


-2.900000 


-3.000000 



TABLE IV: Eigenvalues for j = 3/2, m = ±3/2 and P = -1. 



p = +l 






aw 








R = pujio 


0.1 


0.2 


0.3 


0.4 


0.5 


0.0 


1.974582 


1.951771 


1.931737 


1.914653 


1 


900690 




2.027860 


2.058002 


2.090270 


2.124517 


2 


160601 


0.20 


1.975854 


1.954168 


1.935077 


1.918712 


1 


905206 




2.026477 


2.055157 


2.085916 


2.118632 


2 


153187 


0.40 


1.977225 


1.956967 


1.939328 


1.924407 


1 


912300 




2.025190 


2.052692 


2.082402 


2.114216 




1 A QflQQ 

14oUoo 


0.60 


1.978695 


1.960169 


1.944501 


1.931764 


1 


922023 




z.uzoyyy 


Z.U0U0U4 


o 070790 




2 


145105 


0.80 


1.980265 


1.963777 


1.950603 


1.940800 


1 


934412 




2.022904 


2.048890 


2.077863 


2.109722 


2 


144362 


1.00 


1.981935 


1.967793 


1.957641 


1.951530 


1 


949490 




2.021904 


2.047548 


2.076820 


2.109598 


2 


145751 




0.6 


0.7 


0.8 


0.9 


1 





0.00 


1.890016 


1.882792 


1.879170 


1.879284 


1 


883249 




2.198388 


2.237750 


2.278567 


2.320724 


2 


364112 


0.20 


1.894692 


1.887297 


1.883142 


1.882341 


1 


884998 




2.189469 


2.227369 


2.266780 


2.307602 


2 


349738 


0.40 


1.903101 


1.896895 


1.893767 


1.893790 


1 


897033 




2.183752 


2.221273 


2.260500 


2.301337 


2 


343692 


0.60 


1.915337 


1.911754 


1.911315 


1.914052 


1 


919984 




2.181178 


2.219372 


2.259589 


2.301732 


2 


345703 


0.80 


1.931468 


1.931986 


1.935967 


1.943401 


1 


954263 




2.181673 




2.263861 


2.308514 


2 


355390 


1.00 


1.951530 


1.957641 


1.967793 


1.981935 


2 


000000 




2.185144 


2.227636 


2.273085 


2.321347 


2 


372281 



TABLE V: Eigenvalues for j = 3/2, m = ±1/2 and P = +1. 



p = -1 






auj 






R = ii/ui 


0.1 


0.2 


0.3 


0.4 


0.5 


0.0 


-1.974582 


-1.951771 


-1.931737 


-1.914653 


-1.900690 




-2.027860 


-2.058002 


-2.090270 


-2.124517 


-2.160601 


0.2 


-1.973408 


-1.949771 


-1.929299 


-1.912203 


-1.898688 




-2.029340 


-2.061228 


-2.095471 


-2.131886 


-2.170299 


0.4 


-1.972333 


-1.948166 


-1.927750 


-1.911326 


-1.899125 




-2.030917 


-2.064837 


-2.101522 


-2.140746 


-2.182294 


0.6 


-1.971356 


-1.946953 


-1.927075 


-1.911986 


-1.901920 




-2.032590 


-2.068830 


-2.108428 


-2.151104 


-2.196590 


0.8 


-1.970476 


-1.946126 


-1.927259 


-1.914140 


-1.906981 




-2.034361 


-2.073209 


-2.116190 


-2.162958 


-2.213182 


1.0 


-1.969694 


-1.945683 


-1.928286 


-1.917745 


-1.914214 




-2.036229 


-2.077973 


-2.124808 


-2.176305 


-2.232051 




0.6 


0.7 


0.8 


0.9 


1.0 


0.0 


-1.890016 


-1.882792 


-1.879170 


-1.879284 


-1.883249 




-2.198388 


-2.237750 


-2.278567 


-2.320724 


-2.364112 


0.2 


-1.888950 


-1.883170 


-1.881506 


-1.884089 


-1.891014 




-2.210546 


-2.252476 


-2.295943 


-2.340814 


-2.386963 


0.4 


-1.891353 


-1.888185 


-1.889758 


-1.896163 


-1.907439 




-2.225963 


-2.271567 


-2.318931 


-2.367895 


-2.418308 


0.6 


-1.897068 


-1.897568 


-1.903499 


-1.914877 


-1.931650 




-2.244637 


-2.295011 


-2.347498 


-2.401899 


-2.458034 


0.8 


-1.905924 


-1.911032 


-1.922289 


-1.939595 


-1.962780 




-2.266546 


-2.322758 


-2.381549 


-2.442671 


-2.505900 


1.0 


-1.917745 


-1.928286 


-1.945683 


-1.969694 


-2.000000 




-2.291647 


-2.354724 


-2.420937 


-2.489975 


-2.561553 



TABLE VL Eigenvalues for j = 3/2, m = ±1/2 and P = -1. 



21 



ri 





1 


2 


k 


3 




4 




1 


-1 


3.33287 X 10"' 


1.66774 X 10" 




-1.82477 X 10 




2.69509 X 10" 






1 


-0.333356 


7.8727 X 10" 


5 


4.68693 X 10" 


5 


-7.76405 X 10 


— 5 


2 





2.90019 X 10"" 


7.32430 X 10" 




7.62158 X 10" 




-2.00106 X 10 


— 3 







8.39768 X 10 


0.0737528 




—7.84339 X 10 


-4 


8.33219 X 10 


4 


3 





-3.56377 X IQ-* 


2.97546 X 10" 


-2 


-5.99395 X 10 


-3 


4.35512 X 10" 


3 







-1.55588 X 10"* 


-0.0288713 




8.79491 X 10" 


3 


-2.42401 X 10 


-3 


4 





-1.87177 X 10"^ 


1.61106 X 10" 


3 


4.59520 X 10" 


3 


-8.23406 X 10 


-3 







1.33331 X 10"* 


-8.50242 X 10 


-4 


4.21954 X 10" 


3 


-1.17962 X 10 


-3 



TABLE VII: Fitting coefficients for j = 1/2, m = ±1/2 and P = +1. The upper line in each row is for m = 1/2, and the lower 
line for m = —1/2. 



n 





1 


k 

2 




3 


4 




1 


3.33047 X 10"' 


3.31222 X 10"' 


3.98816 X 10" 


3 


-2.03543 X 10"^ 


2.37639 X 10" 


4 




-3.32884 X 10"' 


-3.34729 X 10"' 


1.82192 X 10" 


3 


-1.36801 X 10"^ 


4.66671 X 10" 


4 


2 


-2.92723 X 10"' 


3.12502 X 10"' 


-1.08338 X 10 


-1 


1.89525 X 10"^ 


-2.68764 X 10 


-3 




-3.00825 X 10"' 


3.09671 X 10"' 


-8.99404 X 10 


-2 


1.03932 X 10"^ 


-3.18256 X 10 


-3 


3 


-8.06550 X 10"^ 


6.85237 X 10"'^ 


9.79197 X 10" 


2 


-6.14738 X 10"^ 


9.59822 X 10" 


3 




8.14284 X 10"^ 


-8.29826 X 10"^ 


3.77385 X 10" 


2 


-2.10873 X 10"^ 


7.55711 X 10" 


3 


4 


4.03322 X 10"^ 


-5.06553 X 10"^ 


-6.96088 X 10 


-3 


2.74214 X 10"^ 


-8.09157 X 10 


-3 




-9.26754 X 10"^ 


4.99333 X 10"^ 


2.46876 X 10" 


3 


1.95452 X 10"^ 


-2.47130 X 10 


-3 



TABLE VIII: Fitting coefficients for j = l/2,m = ±1/2 and P = —1. The upper line in each row is for m = 1/2, and the lower 
line for m = —1/2. 



n 





1 




2 


k 


3 


4 




1 


-1 


1.99945 X 10" 


i 


3.71992 X 10" 


4 


-7.87532 X 10"" 


4.65803 X 10" 


4 




1 


-2.00063 X 10 


-1 


2.84647 X 10" 


4 


-3.91985 X 10"" 


1.89222 X 10" 


4 


2 





2.87850 X 10" 


4 


2.99573 X 10" 


2 


4.57896 X 10"^ 


-2.73991 X 10 


-3 







3.72905 X 10" 


4 


3.02612 X 10" 


2 


2.33862 X 10"^ 


-1.13402 X 10 


-3 


3 





-4.81935 X 10 


-4 


1.25710 X 10" 


2 


-8.94371 X 10"-^ 


4.94845 X 10" 


3 







-6.35154 X 10 


-4 


-6.12322 X 10 


-3 


-3.03913 X 10"^ 


1.83562 X 10" 


3 


4 





2.62158 X 10" 


4 


-1.84259 X 10 


-3 


7.31973 X 10"^ 


-3.94239 X 10 


-3 







3.25972 X 10" 


4 


-1.57711 X 10 


-3 


3.59247 X 10"^ 


-1.49975 X 10 


-3 



TABLE IX: Fitting coefficients for j = 3/2, m = ±3/2 and P = +1. The upper line in each row is for m = 3/2, and the lower 
line for m = —3/2. 



n 





1 




k 

2 




3 


4 


1 


5.99797 X 10"' 


2.00423 X 10" 


1 


3.40583 X 10" 


5 


-4.77152 X 10"" 


2.28612 X 10"" 




-6.00024 X 10"' 


-2.00033 X 10 


-1 


2.87207 X 10" 


5 


7.14085 X 10"^ 


-6.24219 X 10"^ 


2 


-1.26255 X 10"' 


1.24048 X 10" 


1 


-3.12219 X 10 


-2 


2.82263 X 10"^ 


-1.47819 X 10"^ 




-1.27827 X 10"' 


1.28410 X 10" 


1 


-3.22650 X 10 


-2 


-8.40480 X 10"" 


6.80416 X 10"" 


3 


-3.53405 X 10"^ 


3.91967 X 10" 


2 


-9.37469 X 10 


-3 


-5.64216 X 10"^ 


3.06728 X 10"^ 




3.03746 X 10"^ 


-2.96270 X 10 


-2 


6.62180 X 10" 


3 


2.07091 X 10"^ 


-1.47298 X 10"^ 


4 


2.41821 X 10"" 


-1.29284 X 10 


-2 


1.51332 X 10" 


2 


-3.26087 X 10"^ 


-9.82596 X 10"" 




-4.08651 X 10"^ 


3.00329 X 10" 


3 


2.36225 X 10" 


3 


-3.51752 X 10"^ 


1.39410 X 10"^ 



TABLE X: Fitting coefficients for j = 3/2, m = ±3/2 and P = —1. The upper line in each row is for m = 3/2, and the lower 
line for m = —3/2. 



22 



n 





1 


k 
2 


3 


4 




1 


-3.33047 X 10"' 


6.87003 X 10"^ 


-3.52549 X 10"^ 


1.25094 X 10"^ 


1.90819 X 10" 


4 




3.32884 X 10~^ 


-6.52992 X 10"^ 


-1.56447 X 10~^ 


8.64220 X 10"" 


-2.19611 X 10 


— 4 


2 


2.92722 X 10"^ 


-3.12089 X 10"^ 


1.53920 X 10"^ 


-1.50603 X 10"^ 


5.51959 X 10" 


4 




3.00824 X 10~^ 


—3.09458 X 10~^ 


1.36164 X 10~^ 


—6.81511 X 10~^ 


1.38445 X 10~ 


3 


3 


8.06564 X 10"^ 


-7.56215 X 10"^ 


-8.92653 X 10"^ 


5.52054 X 10"^ 


-6.37204 X 10 


-3 




-8.14264 X 10"^ 


8.25932 X 10~^ 


-3.92363 X 10"^ 


1.47876 X 10"^ 


-4.17156 X 10 


-3 


4 


-4.03332 X 10"^ 


5.10444 X 10"^ 


4.71906 X 10"^ 


-2.46746 X 10~^ 


6.22122 X 10" 


3 




9.26669 X 10"^ 


-4.80496 X 10"^ 


-3.88778 X 10"^ 


2.34438 X 10"^ 


-1.17466 X 10 


-4 



TABLE XI: Fitting coefficients for j = 3/2, m = ±1/2 and P = +1. The upper line in each row is for m = 1/2, and the lower 
line for m = —1/2. 



n 





1 


k 
2 




3 




4 




1 


1.99974 X 10"' 


6.65908 X lO"'' 


2.08259 X 10" 


4 


-2.82760 X 10 


-4 


-5.93627 X 10 


-5 




-1.99866 X 10"^ 


-6.69306 X 10"^ 


9.28649 X 10" 


5 


1.25203 X 10" 


4 


-8.67219 X 10 


-5 


2 


-1.91707 X 10"^ 


1.92687 X 10"^ 


-1.24489 X 10 


-1 


3.25586 X 10" 


3 


2.08439 X 10" 


4 




-1.93239 X 10"^ 


1.94559 X 10"^ 


-1.23591 X 10 


-1 


-8.37842 X 10 


-5 


2.56220 X 10" 


4 


3 


-1.65338 X 10"^ 


1.26573 X 10"^ 


-2.56257 X 10 


-2 


-2.31309 X 10 


-3 


-8.48344 X 10 


-4 




1.90423 X 10"^ 


-2.18047 X 10"^ 


3.85236 X 10" 


2 


-9.44174 X 10 


-3 


1.13187 X 10" 


3 


4 


8.26685 X 10"^ 


-1.46583 X 10"^ 


7.67572 X 10" 


3 


-5.46338 X 10 


-3 


7.20356 X 10" 


3 




1.78362 X 10"^ 


-4.28852 X 10"^ 


2.04973 X 10" 


3 


-3.65754 X 10 


-3 


1.31285 X 10" 


3 



TABLE XII: Fitting coefficients for j = 3/2, m = ±1/2 and P = —1. The upper line in each row is for m = 1/2, and the lower 
line for m = —1/2. 



